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The Nambu - Bogoliubov - de Gennes method is applied to the problem of superconducting QCD. 
The effective quark-quark interaction is described within the framework of the Nambu - Jona-Lasinio 
model. The details of the phase diagram are given as a function of the strength of the quark-quark 
coupling constant G' . It is find that there is no superconducting phase transition when one uses the 
relation between the coupling constants G ' and G of the Nambu - Jona-Lasinio model which follows 
from the Fierz transformation. However, for other values of G' one can find a rich phase structure 
containing both the chiral and the superconducting phase transitions. 



CO . I- INTRODUCTION 

o 
o 

In this letter we apply the Nambu - Bogoliubov - de Gennes (NBG) [1] method to the problem of superconducting 
QCD [2,3]. The Nambu - Jona-Lasinio hamiltonian for two flavors and three colors [4] is chosen as a model for the 
I ^ 1 quark interaction: 
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where m is a current quark mass, f is the vector of Pauli matrices, and all indices are suppressed. The hamiltonian 
' (1) describes the quark interaction. It works under the assumption that the deconfinement phase transition appears 
, at a lower baryon density then the chiral phase transition. In that case we can think of quarks as relevant degrees of 
7— I ' freedom at appropriate densities. 

The NBG method is a very convenient tool for mean-field calculations of any hamiltonian containing four-fermion 
interactions. It allows us to trace the physical details of the considered problem and is easy to apply. The problem 
of the superconducting QCD within the framework of Nambu - Jona-Lasinio model has been considered before in [5] . 
This letter differs in the approach and in details of the regularization procedure. Instead of using the three-momentum 
cut-off for all interactions we choose the form-factor regularization for diquark interactions and leave the cut-off only 
to regularize the infinite contribution from the Dirac sea of antiquarks. The cut-off regularization influences only the 
features of the chiral phase transition and has only minor meaning for the superconducting phase transition. 

In sections II and III we describe the Fierz transformation and the mean-field approximation. Section IV gives the 
• details of NBG approach. The last section contains numerical results. 



II. FIERZ TRANSFORMATION 

Through the use of appropriate Fierz transformation one can rearrange the hamiltonian (1) to the form which 
contains color 6 and 3 operators explicitly. To this end, one can use the relations between Pauli matrices: 

TikVj + 5 ik Slj = -{T2)il(T 2 )jk + (fT 2 ) l l(T 2 T) J k 

TtkVj — SikSlj = —2(T 2 )il(T 2 ) : jk (2) 

and the following decomposition of unity in the SU (3) color group: 

5 a f}5 1 s = -5 ai 8 sp + 2 J J . (3) 

From the Fierz transformed form of (1) one can verify that the interactions in the scalar or axial- vector 3 and tensor 
or vector 6 channels are attractive, whereas in pseudoscalar or vector 3 and axial-vector 6 channels are repulsive. 
Assuming that only scalar 3 creates condensate one obtains the relation G' — G/4 between the coupling constant G" 
of the diquark scalar channel and the coupling constant G of the Nambu - Jona-Lasinio hamiltonian (1). Nevertheless 
in a more general treatment one can try to use G" as an additional parameter, as we do in this letter. Thus our 
starting point is the following hamiltonian operator: 
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H = J d 3 x [v£ (-H ■ V + m) V£ - G (& a ipi) 2 - 

G' (lPi(T2) jk (t A ) af3 C l5 4) 1 ^T(T2)mn(t A )s p C l5 r p )] , (4) 

where j, k,m,n = u,d are flavor, a, {3,5, p = R, G, B are color indices and spinor indices are suppressed. The quantity 
t 2 is antisymmetric Pauli matrix, t A is vector of three antisymmetric Gell-Mann matrices and C = 17270 is charge 
conjugation matrix. For G, G' > hamiltonian (4) describes the attraction of quarks in scalar color singlet and scalar 
color antitriplet channels. 

III. MEAN-FIELD ANALYSIS 

To study the finite density effects at zero temperature we have to consider the effective hamiltonian operator, 
F = H — fj,N, where [i is teh quark chemical potential and N is particle number operator: 

N = J d 3 x(^J^ a . (5) 

The non-zero density of quarks influences the structure of the vacuum by developing the energy eigenstates ladder of 
the Fermi sea quarks. However, as is well known, the Fermi sea is unstable in the presence of the attractive channels 
leading to the Cooper instability. Thus we expect the creation of uniform condensates in color singlet and antitriplet 
channels. The non-zero color singlet condensate, called the chiral condensate, breaks the chiral symmetry and leads to 
the phenomenon of mass generation for quarks [4]. The non-zero color antitriplet condensate breaks the color gauge 
symmetry and leads to the phenomenon of Colour Superconductivity [2,3]. Both phenomena can be described in the 
mean-field approximation. Assuming that the fluctuations of the fields around the uniform condensates are small we 
linearize the four-fermion operators: 

(^T 2 t A C^ (^T 2 t A C 75 ^) W (i>T 2 t A C l5 ^ (i>T 2 t A C l5 lP) + 

((^T 2 t A C l5 ^)ijT 2 t A C l5 ^ ~ {(i>T 2 t A C lh ^){i>T 2 t A C l5 i>), (6) 

and 

(#) 2 «2(#)#-{#) 2 . (7) 

The brackets in the above formulae denotes the vacuum expectation values of the field operators. Assuming that the 
lowest energy state is uniform 1 we introduce the chiral condensate M: 



and the superconducting condensate A: 
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where the direction in the color space is chosen arbitrarily. With the new variables, the effective hamiltonian in the 
mean field approximation takes the form: 



F = J d 3 x ^ (-iS ■ V + M 7o - 1> ~ j (ipT^Cjsi)) 1 - ^-ipT^C^ + ^7 + 



(M - to) 2 



AG 

where M — M + m is the constituent quark mass. 



(10) 



1 This assumption is not obvious, see e.g. [8] 
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IV. NAMBU-BOGOLIUBOV-DE GENNES APPROACH 



We use the NBG approach which is particularly useful in practice and enlights the physical aspects of the considered 
problem. To start with, it is very convenient to introduce the creation and annihilation operators: 

/<i 3 fc r - - ~\ 

(2^)3 V2£^ .(fc)expHfc*) +v s {k)V a ] s {k)eMikx)\ , (11) 

where Ek — + M 2 , u s ,v s are Dirac bispinors, a % a s (k) (b^ s (k)) is an annihilation operator of quark (antiquark) 
of color a, flavor j, spin s, and momentum k, satisfying anticommutation relations: 

{< s (k),a^ r (p}} = {b^ s {k)^] r (p}) = (2^f8 sr 5 aP 8(k-p). (12) 

The theories with point four-fermion interactions are not renormalizable and require regularization which introduces 
additional parameters. We use the form-factor regularization, as well as the standard three-momentum cut-off. The 
form-factor regularization introduces the function f(p) for each fermi field in (10) that cuts the large momenta in the 
interaction. Following [3] we use: 

M = ^ (13) 

where oj is of order of 0.4-1.0 GcV. The use of the form-factor is a convenient regularization of color superconducting 
part of (10). However, in the description of chiral symmetry braking it unnecessarily complicates the shape of the Fermi 
surface. Thus for that part we use the three-momentum cut-off regularization. This introduces another parameter, 
A. The parameters G = 5.01 GeV~ 2 , A = 0.65 GeV are fixed within the Nambu - Jona-Lasinio model. They directly 
influence the features of the chiral symmetry breaking [6]. In the numerical calculations we put arbitrarily u> = 0.8 
GeV. The features of the chiral symmetry breaking are robust against the details of the form-factor regularization. 
However, such quantities as the value of the gap parameter A, or the exact position of the superconducting phase 
transition in the phase diagram do depend on these details. We discuss this point in the last chapter. By using the 
decomposition (11) we can write the effective hamiltonian (10) in the form 2 : 



f d 3 p J 
J (2^) 3 \ 
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A.f = Af(p) 2 (r 2 ) jk (t 1 U- (15) 



The symbols |, J. denote the spin projection of (anti)quarks on the direction p. From the formula (14) it follows that 
at the mean field level the quarks and antiquarks decouple from each other and can be treated independently. Thus 
in the following we consider only the quark part of the effective hamiltonian. The first term and the second line of 
(14) describing the interaction of quarks can be write in the form: 



WlH<>-*.<-»)( v ) +2 - s /w ( *-"'- (16) 



2 Time play no role in our considerations thus was chosen arbitrary to t — 0. 
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where the last integral cancels the infinite contribution from the inverse ordering of creation and annihilation operators 
in (16) with respect to (14). The degeneracy factor 2 • 3 counts flavors and colors. The task is to diagonalize the 
hamiltonian F part . It can be checked that (16) can be written down as the sum over six independent families: 
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The symbols R, G, B denote colors and u, d flavors quantum numbers of quarks. The last two families describe the 
"Blue" quarks which are free and do not contribute to the color superconducting phase. The diagonalization of F part 
is now straightforward because the problem has been reduced to the independent diagonalization of 2 by 2 hermitian 
matrices A and B. To this end one introduces the Bogoliubov unitary transformation of the vectors (19). For families 
a = 1, 2 this transformation takes the form: 



U = 



2T ex P(^) 



-W4j^exp(-iJ) rfjsg 



, A = ^+W7W , A = |A|exp(t<y), 



(20) 



where 77,77' are arbitrary phase factors. For families a — 3,4 Bogoliubov transformation is the same as (20) but the 
change of the sign A — ► —A. Through the transformation (20) one can introduce the quasiparticles basis: 
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In the limit of vanishing interaction |A| 
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(21) 



(22) 



The above equations confirm the interpretation of a\ , /?j as the creation operators of particles above the Fermi surface 
and of holes below the Fermi surface. The hole can be understood as a lack of particle from the occupied Fermi sea. 
In the quasiparticles basis F part is diagonal and takes the form: 



F part = [ £ ^e 2 + |A| V(P) 4 K(PK(P) + f3l(p)Pa(P)) 

J ^ Z7T > a=l 



4^ ^ s ea Bjp)<jp) 



4/|^(e-V?T^W). 



(23) 



The factor 4 in the last term comes from 2 flavors and 2 paired colors (except for the "Blue" one). 

The contribution from antiquarks can be find in exactly the same way. The only difference is the Dirac sea 
contribution which is crucial for the chiral symmetry restoration (e.g. [7]). Indeed V a operators in the first line of (14) 
have to be normally ordered. This ordering leave the infinite Dirac sea contribution to the effective hamiltonian. The 
final expression for the full diagonal free energy operator is: 
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(24) 



The 1st and 2nd line of (24) describe the spectrum of quasiparticles constructed from the quarks (a a ,f3 a ) and an- 
tiquaries (a a ,/3 a ). An example of such a spectrum is pictured in the Fig. 1 as a function of \p\. The lower branch 
describes the "quark" and the higher "antiquark" quasiparticle dispersion relations. The minimum of the "quark" 
spectrum coincides with the Fermi energy (fi = 0.38 GeV in the Fig. 1) and its value is equall to the size of the 
energy gap. The low energy exictations around the Fermi surface thus consist of the "quark" quasiparticles and the 
contribution from antiquarks can be neglected. This is the main reason why in the first approximation the relativistic 
theory behaves very much like the non-relativistic field theory of condensed matter system. Due to the Pauli principle 
the pair creations are impotrant only at the energy scales of the order of the double chemical potential /i or in the 
subtle effects like the decay of the " pion" [9] . The 3rd line of (24) describes the contribution from unpaired " Blue" 
quarks and antiquarks. The operator s is defined by the following "Bogoliubov transformation": 



J 



B,s 



E p < (i 
E p > /! 



(25) 



The operator part of the hamiltonian (24) is positively defined. This let us to define the vacuum state |3>) as the 
"trivial" with respect to the new basis operators: 



«„!*> = /?«!*> = 5a|*> = &|*> 



*4» = o 



(26) 



The c-number part of the hamiltonian (24) gives the vacuum energy. In the 5th line we have the contribution from 
pairing interactions of quarks and the last line describes the contribution from the Fermi sea of "Blue" quarks and 
from the Dirac sea (the last term) of all antiquarks. The Dirac sea term is only one divergent in (24). In that place we 
used additional cut-off rcgularization which was discussed before. The constant C is choosen such that the preassurc 
and number of particles at zero density is zero. The gap equations for M and | A| is obtained by the variation of the 
vacuum expectation value of (24) with respect to M and A. 
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FIG. 1. Spectrum of quasiparticles at jx — 0.38 GeV for massive quarks m — 7 MeV and G = G' — 5.01 GeV -2 
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V. RESULTS 



Let us first consider the case of massless quarks m — 0. The values of G = 5.01 GcV~ 2 and A = 0.65 GeV can be 
fitted to the values of the chiral condensate and the pion decay constant [6]. This choice of parameters establishes 
that the chiral phase transition is first order, its position is at fi = 0.316 GeV, and the size of constituent quark mass 
at the transition point is M = 0.29 GeV. The values of G', u have little influence on these results. On the other hand, 
the position of the superconducting phase transition (CS) and the size of diquark condensate A do depend on these 
parameters. The value of lo — 0.8 GeV has been chosen arbitrarily and the behavior of CS phase transition has been 
checked as a function of G'. For low values (numerically below G/2) the chiral and CS phases are separated in the 
density phase diagram. At higher values of G' there are coexistence regions. Let us consider two cases G' = G/4 
and G' = G. For G' = G/4 the chiral phase transition is placed at /j, = 0.316 GeV and then at /j, = 0.7 GeV only 
a trace of continuous CS phase transition appears (in a sense that the value of A reach roughly 1 MeV at best and 
quickly vanishes). The vanishing of A at high density is a general feature of the model and is discussed at the end 
of this section. For G' — G another scenario is realized. First at around fi = 0.305 GeV there appears a continuous 
phase transition to the CS phase at non-zero M = 0.301 GeV. The gap A develops up to the value of 3 MeV at 
fi w 0.317 GeV where there is the first order chiral phase transition. The gap parameter A changes from about 3 
MeV to 18 MeV. The size of the gap A is not very high in compare to e.g. [3] however it size depends on the values 
of the parameters G',uj. Increasing the density increases the size of A. It reaches a maximum value of A = 43 MeV 
at // = 0.55 GeV. 

In the real world the chiral symmetry is explicitly broken by current masses of quarks. For u, d quarks we take them 
equal, m = m u = nid = 7 MeV. Inclusion of massive quarks leads to a little bit higher values of the gap parameter 
A in the coexistence phase compared to the massless case and changes the position of phase transitions in the phase 
diagram. As an example let us consider the case of G' = G. The phase diagram is given in the Fig. 2. It can be seen 
that at At ~ 0.33 GeV we have a continuous superconducting phase transition. The value of the gap A increases with 
increasing density up to about 8 MeV at At ~ 0.341 GeV. Then there appears a first order chiral phase transition and 
parameters change from M = 0.29 GeV to M = 0.15 GeV and A = 8 MeV to A = 20 MeV (in the Fig. 3 there is 
the contour plot of the potential as a function of M, A at At = 0.341 GeV). When the density is further increased the 
value of M decreases and A increases with maximum around 43 MeV at fi w 0.56 GeV. 

For both massive and massless quark cases at high enough densities (a* around 0.9 GeV) the gap A vanishes. This 
comes as a result of the regularization by the form-factor. At high densities the energy of interaction is too small to 
overcome the volume energy of the creation of the superconducting condensate. In this region the superconducting 
gap has to be calculated through the use of another approximation. 

In conclusion, one can find that within the NJL model, using the relation G' = G/4 between diquark and quark- 
antiquark channels coupling constants, there is no superconducting phase transition in the region of model applicability. 
However, when one releases this constraint, and varies the value of G', there are rich structures of the chiral and 
superconducting phase transitions. For lower values of G', these transitions are separated in density phase diagram, 
whereas for higher values there is a coexistence region. The non-zero current quark masses has only minor influences 
for the whole picture. 
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FIG. 2. Dependence of M and A on chemical potential fi 
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FIG. 3. Contour plot of free energy as a function of M and A at the vicinity of the phase transition at fi = 0.341 GeV. 
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